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With the advent of high analysis technology in detecting the Gamow-Teller (GT) 
excited states beyond one nucleon emission threshold, the quenching of the GT 
strength to the Ikeda sum rule seems to be recovered by the high-lying GT states. 
Moreover, in some nuclei, the stronger GT peaks than any other peaks appear ex- 
plicitly in the high- lying excited states. We address that these high- lying GT excited 
states stems from the smearing of the Fermi surface by the increase of the chemical 
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potential owing to the deformation within a framework of the deformed quasi-particle 
random phase approximation (DQRPA). Detailed mechanism leading to the smear- 
ing is discussed, and comparisons to the available experimental data are shown to 
explain the strong peaks on the high-lying GT excited states in a satisfactory manner. 

PACS numbers: 23.40.Hc, 21.60.Jz, 26.50.+X 



^ ! The deformation in nuclei becomes more important than last decades with the recent 

Sh ■ 

development of rare isotope (RI) accelerator facilities, from which one may perform lots of 
challenging experiments related to the RI nuclei. Most of them are thought to be deformed 
more or less. They are usually produced in the successive nuclear capture reactions in the 
cosmos, i.e. slow- and rapid- process, presumed to be occurred, respectively, at the initial 
and explosive stage of stellar evolution, such as supernovae (SNe) explosion. Although they 
decay fractions of a second, they imprint their existence on the nuclear abundances of stars 
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Another interesting nuclear process is the neutrino (z/)-process in the SNe explosion which 
is treated as one of important nuclear processes for the nucleosynthesis, in particular, of the 
neutron-deficient nuclei, such as 138 La and 180 Ta [2-4 1. Here 180 Ta is a well known deformed 
nucleus. Since the neutrinos emitted from a proto-neutron star may have tens of MeV energy 
high enough to excite the deformed nuclei, one needs to understand more precisely the high- 
lying excited states beyond one nucleon threshold. Among them, the GT excitations are 
of great importance because most of the charge exchange reactions are dominated by the 
GT transition. These high-lying GT excited states are also closely related to the nuclear 
structure. For example, the GT quenching problem says that the difference of total running 
sums S± for empirical GT (±) transitions, 5_ — S + , is usually quenched compared to the 
Ikeda sum rule (ISR), S_ — S + = 3(N — Z). 

However, recent experimental GT data on the high-lying states deduced by the more 
energetic projectiles shed a new light on the GT states located above one nucleon threshold, 
whose contributions are thought to enable us to explain the quenching problem with the 
ligh-lying GT excited states through the multi-particle and multi-hole configuration mixing 
5)]. Of course, the contributions from the A excitation and from the two-body current may 
be still effective for the quenching problem. 

Conventional approach to understand the nuclear structure is based on the spherical 
symmetry. In order to describe neutron-rich nuclei and their relevant nuclear reactions 
occurred in the nuclear processes, one needs to develop theoretical formalism including 
explicitly the deformation Ref. Q exploited the Nilsson basis to the deformed quasi- 

particle random phase approximation (DQRPA). But the two-body interaction was derived 
from the effective separable force. The realistic two-body interaction is firstly considered at 
Ref. [3] only with neutron-neutron (nn) and proton-proton (pp) pairing correlations which 
have only isospin T = 1 and J = interaction. But, to properly describe the deformed 
nuclei, the T= and J = 1 pairing should be also taken into account because it may easily 
lead to the deformation by the J = 1 pairing. 

In this work, we extend our previous QRPA based on the spherical symmetry jsj], which 
has been exploited as a useful framework for describing the nuclear reactions sensitive on 
the nuclear structure of medium-heavy and heavy nuclei [9]. Since our DQRPA formalism 
is fully discussed at Ref. [l0|, we briefly summarize the DQRPA used in this work. Starting 
from the deformed Wood Saxon potential we transform a physical state given by the 
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diagonalization of total Hamiltonian in the Nilsson basis into the spherical basis, in which 
one may perform more easily theoretical calculations. 

In a cylindrical coordinate, eigenfunctions of a single particle state and its time-reversed 
state denoted as a and a in deformed Woods-Saxon potential are expressed as follows 



\a Pa = +1 >= Y^N^Oa \ N > n *> A «' fi « = A « + 1/2 > 

Nn z 

+&kinJW,n*,A a +l,fl Q = A a + 1- 1/2 >], (1) 



where A" = nj_ + n z (n± = 2n p + A a ) is a major shell number, and n z and n p are numbers 
of nodes of the physical state on the deformed harmonic oscillator wave functions in z and 
p directions, respectively. A a (Q a ) is the projection of the orbital (total) angular momentum 
onto the nuclear symmetric axis z. The coefficients q are obtained by the eigenvalue 
equation of the total Hamiltonian in the Nilsson basis. 

The deformed harmonic oscillator wave function, \Nn z A a ll (= ±1/2) > in Eq. flTJ can 
be expanded in terms of the spherical harmonic oscillator wave function \NqIAH > as follows 

AT J —Nq—1 

\Nn z A a Z>=J2 Ktr^ £ °£iE W ^ >> ( 2 ) 
AW i 

with the Clebsch-Gordan coefficient C^°i s and the spatial overlap integral Aj^ A =< 
N lA\Nn z A > numerically calculated in the spherical coordinate system. Therefore, the 
expansion can be simply written as \aQ a >= ^ a B% \aQ a >, with the expansion coefficient 
summarized as B« = J2 Nn ^ C^A^ 

The Deformed BCS equation is solved by using ab initio Brueckner G-matrix calculated 
from the realistic Bonn CD potential for the nucleon-nucleon interaction. For example, the 
pairing potentials A p is calculated as 

A - =-- - Vo p .F J0 F J0 

^apap 2 (2j _j_ I) 1 / 2 / 4 Hpair ± aaaa 1 7C7C 

G{aacc, J)(2j c + l) 1/2 {u* lp v lpc + u* 2p v 2pc ) , (3) 

where = B^ Cj% a jga.SK' = Q a + is designed for the transformation to 

the deformed basis of the G-matrix. Here K', which is a projection number of the total 
angular momentum J onto the z axis, is selected K' — at the BCS stage because we 
consider the pairings of the quasi-particles at a and a states. In order to renormalize the G- 
matrix, strength parameters, g pair , g pair and g p 2 ir are multiplied to the G-matrix to reproduce 



empirical pairing gaps [8]. The ^ decay operator, fif 1 , is defined in the intrinsic frame as 

= < a pPa\ T+(J K\PnPl3 > C ] p C n . (4) 

The (3 transition amplitudes from the ground state of an initial nucleus to the excited state 
are expressed by 

< l(K),m0-\ QRPA > 

= N aa " Pa pf3"pp < OKx'-pp a \oK\fifi"npp > 

aa"p a W"pp 

[Upaa" v nPP"Xaa"PP",K + v paa"U n pi3"Yaa"PP",K] , (5) 

where | QRPA > denotes the correlated QRPA ground state in the intrinsic frame with 



the nomalization factor given as J\f aa "j3p"(J) = a/1 — Saf3&a"i3"{— 1) J+T /(1 + 5 Q/ 3<5 a ",3")- The 
Wigner functions disappeared by using the orthogonality of two Wigner functions from the 
operator and the excited state, respectively. To compare to the experimental data, the 
GT(±) strength functions 

4H= E \<l(K),m\\P±\\QRPA>\ 2 (6) 

K=0,±1 

and their running sums S^ T = £ m B± T (m) for 76 Ge and 82 Se are enumerated ratio of 
the Sq T — Sq T to the ISR value 3(iV — Z) calculated with the closure relation (ISR I) and 
without it (ISR II). 

The single particle states are used up to 4hu for two nuclei, in the spherical limit. Since 



the GT strength distribution turns out to rely on the deformation parameter, /3 2 [111 ], we 
exploited the deformation parameter, |/3 2 | < 0.35, as the default values. Corresponding 
values of the deformation parameter /3 2 , the renormalized strengths g pair and gt air , and 
theoretical and empirical pairing gaps are tabulated in Table I. The ISR I and II, which are 
deviated 8% maximally from the ISR, seems to be reasonable more or less, if we remind 
that the particle model spaces spanned by the deformed basis are finite less than 4 fuv and 
two-body currents including A-excitation are not taken into account in this work. 

An interesting and important point to be noticed is that the ir values are abnormally 
deviated from 1 for some prolate /3 2 values. This turns out to results from the wide smearing 
of the Fermi surface in deformed nuclei, which makes the uv coefficients for the paring gaps 
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TABLE I: Deformation parameters fa, pairing strength parameters 9^ a i r (gp a i r ) , and theoretical 
pairing gaps A^™ used in this work. The ISR I and II in the last column are given as (Sq T — 
Sq T )/3(N — Z). The particle-particle (particle-hole) strength paramet ers arc taken as Qpp — 
1.0(g ph = 1.0) for two nuclei. Empirical pairing gaps, A p em (A" m )=1.561 (1.535) and 1.409 (1.544) 
Mev, are exploited for 76 Ge and 82 Se, respectively. 
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fa 


<f ■ (q n ■ ) 
n pair \n pair J 




An 


ISR I, II(%) 


76 Ge 


0.1 


1.10(1.37) 


1.575 


1.538 


97.04, 106.03 




0.2 


1.17(1.41) 


1.578 


1.540 


96.94, 102.40 




0.35 


1.52(3.19) 


1.562 


1.545 


96.48, 99.54 




0.1 


1.01(1.49) 


1.502 


1.579 


96.51, 106.45 


82 Se 


0.2 


1.11(3.49) 


1.419 


1.547 


95.62, 107.55 




-0.1 


1.02(1.53) 


1.410 


1.564 


96.58, 104.89 




-0.2 


1.01(1.64) 


1.568 


1.549 


96.57, 109.69 



in Eq.(3) small around the Fermi surface and enlarges the g™ air to fit the empirical pairing 
gap. 

This tendency is explicitly revealed in the occupation probabilities of 82 Se in Fig.l, which 
shows a wide smearing, i. e. significant change of the occupation probabilities of the particles 
around the Fermi surface by the increase of the deformation parameter fa. Since the increase 
of the prolate deformation usually raises the Fermi surface energy, some states above the 
Fermi surface burrowed below or around the Fermi surface. For example, (4 3 1 1/2), (4 2 
2 3/2) and (4 11 3/2) states in 82 Se are such cases. This gives rise to the wide smearing. 
Such wide smearing of the Fermi surface is also found for the protons and neutrons in 76 Ge, 
as shown in Fig. 2. 

Detailed changes of single particle states by the fa deformation leading to the smearing 
are presented in Figs. 3 and 4. Some states above the Fermi surface are reallocated below 
or around the surface, because the single particle state energies adopted from the deformed 
Woods Saxon potential depend on the parameter fa. The deformation of nuclei may be 
conjectured to come from macroscopic phenomena, for example, the core polarization, the 
high spin states and so on. Microscopic reasons may be traced to the tensor force in the 
nucleon-nucleon interaction, which is known to account for the shell evolution according to 
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FIG. 1: (Color online) Occupation probabilities of the neutrons in 82 Se as a function of the single 

particle energy given by Nilsson basis for two deformation parameters /?2- With the increase of /?2 

from 0.1 to 0.2, the Fermi energy A n is increased with the wider smearing of the Fermi surface. 
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FIG. 2: (Color online) Same as Fig.l but with 76 Ge. For comparison, results for protons are shown 
in the left panel. 
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the recent shell model calculations For example, T = 0, J = 1 pairing, which is 

associated with the 3 So tensor force, may lead to the deformation compared to the T = 1, 
J = pairing. 

Therefore, the deformation parameter adopted in this work may include implicitly and 
effectively such effects, because the single particle states from the deformed Wood Saxon 



potential show strong dependence on the 02, as shown in Figs. 3 and 4 14]. We have 
taken g p h(g PP ) = 1.0(1.0) for the particle-hole (particle-particle) strength parameters for 
the two nuclei. Since the nuclei considered here are expected to have large energy gaps 
between proton and neutron spaces, in this work, we considered only the nn and pp pairing 
correlation although the formalism is presented generally. For example, in the neutron-rich 
nuclei of importance in the r-process, the np pairing may not contribute so much. But for the 
p-process, the np pairing could be more important than the neutron-rich nuclei because of 
the adjacent energy gaps of protons and neutrons. The calculations for the neutron-deficient 
nuclei in p-process are in progress. 

In the following, we show the GT strength distributions in terms of the 02 parameter, as 
a function of excited energy of parent nuclei. Therefore, experimental data are presented 
by adding the empirical Q values from the measured data. In Fig 5, the GT strength 
distribution, B(GT~), on 82 Se is presented as a function of the excitation energy E ex w.r.t. 
the ground state of 82 Se for 02 = ~ ±0.2. Uppermost panel represents the experimental 
data deduced from the 82 Se(p,n) 82 Br reaction at 134.4MeV With the redistribution 

of the GT strengths, the GT data around 12MeV are well reproduced at (3 2 = 0.2 (see (c) 
panel in Fig.5). For a reference, the 02 value from the relativistic mean field (RMF) is 
0.133 Q. Fig.6 shows B(GT-) of 76 Ge function of the excitation energy E ex with 
prolate shapes, 0% — ~ 0.35. Uppermost panel represents the experimental data from the 
76 Ge(p,n) 76 As reaction at 134.4MeV 15], which show a strong GT state peak around 12 
MeV. The GT strength distributions are widely scattered due to deformation similarly to 
Fig.5. In particular, results for the 0i = 0.35, in which the ISR I and II are almost satisfied, 
nicely reproduce the peak on the high-lying excited states. The value calculated by RMF is 



0.157 [16|. 

Both nuclei, 76 Ge and 82 Se, show a strong peak on high-lying GT states around 12 MeV. 
Our 02 values, 02=0-35 and 0.2 for 76 Ge and 82 Se, customized to the ISR values, are shown 
to reproduce neatly the GT states. Actually, these high-lying GT states turn out to be 
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FIG. 3: (Color online) Single particle energies of 82 Se as a function of the different deformation 
parameter fa. Red dotted lines denote the Fermi energy. 



intimately associated with the smearing of some physical states around the Fermi surface in 
deformed nuclei. Figs. 7 and 8 show the forward amplitudes X of two excitation energies 
and one excitation energy, having the large B(GT) values in (c) and (d) of Figs. 5 and 6, 
respectively. For 82 Se, the main high-lying GT transition come from (4 3 1 3/2)— » (4 4 
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FIG. 4: (Color online) Same as Fig.3 for 76 Ge 

1/2), (3 1 3/2) ^ (3 1 1/2) and (4 1 1 3/2)-» (4 1 1 1/2) state transitions, while (3 3 
7/2)-K3 3 5/2), (3 1 3/2) (3 1 1/2), (4 1 3 7/2)^ (4 1 3 5/2) and (3 1 2 3/2)^ (3 
1 1/2) turn out to be main transitions for 76 Ge. 

In summary, to describe the single particle state in deformed basis, we used the deformed 
axially symmetric Woods-Saxon potential. We performed the deformed BCS and deformed 
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FIG. 5: (Color online) Gamow- Teller strength distributions B(GT ) on 82 Se as a function of the 
excitation energy E ex w.r.t. the ground state of 82 Se. Experimental data denoted as filled (red) 
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X 

CD 



Q_ 

E 

CO 
"D 
CO 



X 




Q_ 

E 

CO 

TD 
i— 

CO 



(a)E =10.74MeV 

v ' ex 

B(GT)= 2.71 ; 




0.5- 



0.0 



(3013/2, 3011/2) 



(b)E =11.51MeV 
B(GT)= 2.81 



(4113/2, 4111/2) 







50 



100 



config. 



FIG. 7: Main transitions in the forward amplitudes X for two dominant high- lying GT states {E ex 
= 10.74 and 11.51 MeV) w.r.t. the ground state of 82 Se as a function of the configuration state for 
the panel (c) in Fig. 5. 

QRPA with a realistic two-body interaction calculated by ab initio G-matrix based on Bonn 
potential. Results of the Gamow- Teller strength, B(GT), for 76 Ge and 82 Se show that the 
deformation effect leads to a fragmentation of the GT strength and reproduces the high-lying 
GT excitation deduced by higher energy projectiles, which result from the wide smearing by 
the increase of the Fermi surface energy due to the deformation. This work was supported by 



13 



X 1.0 



CD 

"D 



Q. 

E 

CO 
"O 
03 



0.5- 



0.0 



~ i 1 1 x~ 



— I 1 1 p" 



E =11.86MeV 

ex 



(3037/2, 3035/2) 

(3013/2,3011/2) B(GT)=6.45 



(4137/2, 4135/2) 

(3123/2, 3101/2) 



r " — ^r- 







I ' " 'I' 



H "i — u L i ^-r" r- 



50 



100 



config. 
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